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The Atmospheric Processing Module (APM) is a Mars In-Situ Resource Utilization (ISRU) 
technology designed to demonstrate conversion of the Martian atmosphere into methane for 
use as rocket propellant for an ascent vehicle. The Martian atmosphere mainly consists of 
95% carbon dioxide (COz) and residual argon and nitrogen. APM utilized cryocoolers for CO2 
acquisition from a simulated Martian atmosphere and pressure. The captured CO2 was 
sublimated and pressurized as a feedstock into the Sabatier reactor, which converted CO2 and 
hydrogen to methane and water. The Sabatier reaction occurred over a packed bed reactor 
filled with Ru/AlOs catalyst pellets. During performance testing at high flow rates, failure of 
the Sabatier catalyst occurred, which was detected by unwanted byproducts at elevated 
reaction temperatures. A Sabatier catalyst failure study was initiated after this event. This led 
to an effort to develop a model that could determine acceptable flow rates and predict the 
temperature profile in order to avoid catalyst damage and identify system tolerances. 
Validation of a computational model with experimental results also aids in long duration 
system development for reactor scaling purposes. This paper discusses the experimental 
results of the APM Sabatier tests, as well as a small scale Sabatier reactor that was set up for 
further model validation. The thermal and kinetic modeling approaches are discussed, as well 
as results of the APM catalyst characterization which included x-ray powder diffraction and 
x-ray photoelectron spectroscopy. 


Nomenclature 


= Reactor cross-sectional area, m? 
= Aluminum oxide 
= Atmospheric processing module 
= Particle shape factor, dimensionless 
= Biot number, (=h Dy 2" keff! ), dimensionless 
= Heat capacity of mixed gas stream, J mol! K"! 
= Hydraulic diameter, (=Z, — Z,), m 
= Equivalent spherical pellet diameter, m 
= Activation energy, J mol! 
Darcy friction factor, (= (0.79In(Re)-1.64)7), dimensionless 
Gas chromatography 
= Convective heat transfer coefficient, Wm? K"! 
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AHR = Heat of reaction, J mol! 


ID = Inner diameter 
ISRU = In-situ resource utilization 
k = Thermal conductivity of fluid, W m?K! 
Ktea = Stagnant bed thermal conductivity, Wm? K"! 
Keat = Thermal conductivity of catalyst, Wm? K"! 
k. = Equilibrium rate, as a function of temperature, k.(T) = exp[(1.0/1.987)(56,000T? + 34,633T"! 
- 16.4In(T) + 0.00557T) + 33.165]!, atm 
Kerr = Effective radial thermal conductivity, W m7? K"! 
ky = Rate constant, = 491,388.9!, s! atm?!* 
MW = Molecular weight of mixture, kg mol"! 
n = Experimental correction factor from Lundel, = 0.225 
N = Molar flow rate, mol s! 
Nu = Nusselt number, dimensionless 
OD = outer diameter 
p = Partial pressure, atm 
P = Pressure, atm 
Pe = Peclet number, ( = UaygD,pCp/k), dimensionless 
Pr = Prandtl number, ( = Cp u/k), dimensionless 
r = Reaction rate of disappearance of CO2, mol m?s"! 
: = Radial coordinate 
R = Ideal gas constant, J mol! K"! or Pa m3 mol! K! 
Re = Reynolds number, ( = pVaygDp/p and pVaygD,/(1 — €)p), dimensionless 
Ru = Ruthenium 
sccm = Standard cubic centimeters per minute 
T = Temperature, K or °C 
TC = Thermocouple 
Uavg = Average velocity, ms! 
Vere = Effective overall heat transfer coefficient, Wm? K”! 
Vi = Stoichiometric coefficient for species 1, Vy2= -4, Veo2 = -1, Vous = +1, Vu20 = +2 
Vv = Reactor volume variable, m? 
VR = Total Reactor volume, m* 
Xx = Conversion 
XPS = X-ray photoelectron spectroscopy 
XRD = X-ray powder diffraction 
y = mole fraction 
Zy = Inner reaction chamber radius, m 
Z = Outer reaction chamber radius, m 
Greek Letters: 
B = Ratio of the surface area at heat transfer interface to reactor volume, m’! 
€ = Fixed-bed porosity or void fraction 
Xr Wavenumber, nm 
y Effectiveness factor 
0 Mass density, mol m? 
Mm = Dynamic viscosity, mol m!s! 
g = Fluid-to-solid thermal conductivity ratio, dimensionless 
Indices: 
i = species 
0 = initial 
f = feed into reactor 
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I. Introduction 


HE success of sustainable human 

exploration and habitat development 
will rely on the explorers to survive from 
resources that are available from the 
explored planet or surrounding region. This 
concept of ‘living off the land’ is called In- 
Situ Resource Utilization (ISRU). On Mars, 
ISRU can include concepts such as capturing Sabatier 
water that is trapped within regolith, Reactor 
utilizing radiation protection from lava- 
tubes, and using the sunlight to convert solar 
energy into electrical power.”~ In this paper, 
the particular concept of harvesting water 
and fuel from the Mars atmosphere via 
cryocoolers and the Sabatier reaction is 
discussed, with a specific investigation of 
the thermal and kinetic behavior of the 
Sabatier reactor. At Kennedy Space Center (KSC), the Atmospheric Processing Module (APM) first collects CO2 from 
a simulated Martian atmosphere (95.4% CO2, 3.0% No, and 1.6% Ar) and Martian pressure (~8-10 Torr) using dual 
cryocoolers.* The CO> is collected by flowing the Mars simulated atmospheric gas into the cryocooler chamber where 
a copper cold head is cooled to 150 K. The COz that is isolated on the cryocooler cold head at 150 K as solid ice is 
then sublimated and stored at a higher pressure (6.8 atm) before it is introduced into the Sabatier subsystem (~ 2.72 
atm). A 3-dimensional model of the main APM components was made with Parametric Technology Corporation Creo 
Elements software and is displayed in Figure 1. The Sabatier reactor is shown in the enlarged box to the left, and the 
cryocoolers are in the enlarged box to the right. The cryocooler is shown with and without its cover to display the high 
surface area copper cold head fitting. Details of the APM cryocooler and Sabatier subsystem are also discussed in 
detail elsewhere.*° The Sabatier reaction, listed in Equation Eq. (1), is an exothermic reaction (AH = -165.4 kJ/mol) 
where COz is reduced by H2 at elevated temperatures over a selective catalyst to produce CH, and H20. 


APM System 


Cryocoolers 
without and with cover 


Figure 1. APM system model. Left: Sabatier reactor; Right: Cryocoolers. 


CO) + 4H2  CHy + 2H20 (1) 


On APM, the H20 vapor is condensed out of the product stream, collected in a vessel, and weighed. In the future, 
it will be sent to the Water Cleanup Module where it will be cleaned for preparation of water electrolysis, which will 
yield Oz that can be stored and liquefied into cryogenic fuel as well as H2 which can be recycled back to the Sabatier 
subsystem to produce CH4. The CH, can be purified and stored before liquefaction for use as fuel, since the primary 
fuel for a Mars ascent vehicle is LO2 and LCHg.’ The concept of the ISRU propellant production system to fuel a Mars 
ascent vehicle is also discussed in detail elsewhere.* 

Sabatier reactors emit heat and are sensitive to many thermal and system conditions, so thermal control is 
imperative for catalytic efficiency and reaction success. On Mars, the low atmospheric pressure and large temperature 
swings reduce thermal conductance and cause temperature inconsistencies, and thermal controls will be required. 
Autonomous operation of a chemical system such as APM is needed to operate for approximately 500 days to generate 
the required fuel for a return vehicle. The fuel production for an ascent vehicle will require approximately 8-10 times 
more production than the current scale of the APM Sabatier reactor. It is necessary to determine the correct reactor 
volume that produces enough resources, while not overheating or generating unwanted byproducts. Modeling has 
begun to characterize the current Sabatier reactor and compare computational results with experimental data, in order 
to understand how to better design future versions of scaled up systems. Modeling the system will also aid in 
continuous control, as thermal runaway can cause catalyst and system hardware damage. 


II. APM Sabatier Reactor Subsystem and Characterization 


A. The APM Sabatier Subsystem 

A simplified APM fluid schematic of the Sabatier subsystem is displayed in Figure 2. H2 was first fed into the 
reactor at 1000 sccm during a preheat treatment from a room temperature k-bottle, until 210 °C was achieved inside 
of the catalyst bed. This Hz reduction step was necessary for catalyst performance. At 210 °C, Hz flow was increased 
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Figure 2. APM Sabatier subsystem schematic. 
to 3375 sccm and CO? was introduced at 750 sccm, and the subsequent temperature rise indicated that the exothermic 
Sabatier reaction had begun. The flow, temperature, and pressure were observed and controlled with LabVIEW 
software. The initial feed ratio of H2 to CO2 was 4.5:1 to allow excess Hz to recycle throughout the system. After the 
reaction was initiated and a catalyst temperature between 400 and 420 °C was achieved, the H2:CO2 molar ratio was 
reduced to 4:1. The CH and H2O vapor products were separated via a shell and tube heat exchanger (chiller/condenser 
system). Water production was measured by periodically removing water from the water collection vessel and 
measuring the mass. The CH, was separated from excess Hz and unreacted COz with a membrane separator. CH4 
purity was detected with a Varian Micro Gas Chromatograph (GC). The permeate stream containing H2, unreacted 
CO? and a small amount of CH, was recycled back to the Sabatier reactor via a pump to complete the reaction. The 
membrane separator relied on a 1.92 atm pressure differential to produce >99.9% pure CH4. The retentate side of the 
membrane was 2.38 to 2.72 atm and the permeate side was approximately 0.75 atm. This pressure differential was set 
manually with a back pressure regulator and required adjustment during any off nominal reading, but can be automated 
in the future. Unwanted byproducts of the reaction included carbon monoxide (CO), and could be detected by the GC. 
CO forms due to changes in feed ratio conditions or changes in temperature. With nominal temperature and pressure 
conditions, the reactor operated at steady state production with the H2 recycle loop for 7 continuous hours. 

Figure 3 shows a model of the Sabatier Reactor. The reactor is a stainless steel tube that is 30 cm long with an 
outer diameter (OD) of 2.54 cm and a wall thickness of 0.2 cm. The reactor is filled with 88 grams of 0.5 weight 
percent ruthenium supported on aluminum oxide (Ru/A12O3) cylindrical catalyst pellets (3 mm) supplied by Sigma- 
Aldrich. Ru/A1O3 catalyst is a proven catalyst with high CHy selectivity.!® The reactant gas entered a feed pipe in the 
center of the reactor so that some of the excess heat of the reaction was used to preheat the incoming gas. The preheat 
loop was used as a heat sink to remove some of the heat from the reaction and avoid reactor thermal runaway. Four 
protrusions on the main portion of the reactor are fixtures for thermocouples (TC) to measure the temperature inside 
the catalyst bed (TC 6, 8, 9, 10). TC 5 is located on the OD of the reactor tube, and TC 14 was inserted inside the 
preheat tube at varying locations (A, B, C) to determine the preheat tube gas temperature. TC 14 was used only for 
preheat temperature 
measurements, and was Joutue 
removed during 
nominal operation. TC 
15 and TC 7 measured rise 
the inlet and outlet of the = 
preheat tube, 
respectively, and were 
also added only for 
experimental validation — farce 

TC14—A =42cm 


and removed during Figure 3. Cross-sectional view of Sabatier reactor model. 


4 


Catalyst Bed 


y) TC15 


TC14 


TC14—C=20cm 


International Conference on Environmental Systems 


nominal testing so as not to interfere with gas flow. The reactor was wrapped with a layer of heat tape (for preheating) 
and was covered with a layer of insulation to maintain as close to an adiabatic condition as possible. 


B. APM Catalyst Characterization 

In 2015, after several months of successful 
nominal testing, the reactant flow rates were 
purposely varied from nominal conditions in order 
to characterize reactor flow performance limits. 
CO, an unwanted byproduct, was detected by the 
GC in the product stream during a high flow rate 
experiment (1250 sccm COz and 5000 sccm Hz). 
The maximum observed temperature was 586 °C 
which was well above the 400 °C nominal 


temperature value, but under the 600 °C safety | | B. i 3 io a 
shutdown temperature value. Although 600 °C was New Pellet Pellet from Failed Inlet Middle Outlet 
Reactor Bed | Pellets from Used Catalyst Bed —| 


not observed by a thermocouple, thermal hotspots 
higher than 586 °C could have been present. The 
Sabatier reaction was unable to run at nominal 
conditions after this run, so it was determined that the catalyst was damaged. The catalyst was removed (Figure 4) and 
some pellets were broken along with residual powders, and the catalyst had turned gray. The catalyst was replaced in 
the reactor bed with a fresh, unused batch of catalyst pellets and the reactions performed nominally once again. 

A Sabatier catalyst failure study was initiated after this event, and led to the computational modeling development 
effort, with the aim of developing a model that could determine acceptable flow rates and predict temperature profiles 
in order to avoid damaging the catalyst and troubleshoot system tolerances. Five Ru/A1203 pellets were analyzed from 
the APM project during the catalyst investigation. They are shown in the bottom of Figure 4 and included a new pellet 
that was unused, a pellet from the Sabatier reactor that failed in 2015 after overheating and running the reactor at 
higher experimental flow rates, and samples from the current nominal Sabatier reactor in its three temperature zone 
regions: catalyst bed inlet, catalyst bed middle, and catalyst bed outlet. 

X-ray powder diffraction (XRD) patterns were obtained from fresh and used catalysts in air at room temperature, 
using a PANalytical Empyrean machine with Cu Ka radiation (4 = 0.15406 nm) at a scan rate (20) of 0.02° s"!. The 
accelerating voltage and applied current were 40 kilovolts and 40 milliamps, respectively. The samples were first 
prepared by grinding the pellets into powder form. Then, additional pellets from the used (but not damaged) catalyst 
bed inlet, middle, and outlet were prepared by scraping off the black/gray Ru layer and separating it from the white 
Al203 support. All the XRD prepared powders were ground in acetone with a silicon carbide mortar and pestle. The 
solids were extracted from the ground solution with a pipette and placed on the powder sample holder to accumulate 
in an even powder layer and dry in air before analysis. The International Centre for Diffraction Data library was used 
from the PANalytical XRD HighScore software for reference peak validation. 

X-ray photoelectron spectrometry (XPS) measurements were conducted with a Thermo Scientific K-Alpha 
instrument using Al Ka radiation on the whole pellets (not the powders of Ru separated from the Al2O3 support). 
Surface survey scans and depth profiling were collected on the pellets. The depth profiling used a 3000 eV etch energy, 
5 second etch cycle, at ~0.35 nm of etching per second. The reference Cls for this work was 284.91 eV. 


Figure 4. Top: Left: The 2015/spent catalyst that was damaged; Right: 
New/unused catalyst; Bottom: Ru/AI,O; pellets that were characterized. 


III. Small Scale Sabatier Reactor 


Thermal and chemical computational modeling on the APM system with a preheat loop posed more challenges than 
the model of a simple packed bed reactor with no internal preheat tube. A small scale Sabatier reactor was set up and 
experiments were run to aid in model validation and compare with some of the literature results found on packed bed 
reactors. The reactor described earlier (Figure 5) used in APM experiments was denoted as “APM Sabatier reactor” 
in this work while the smaller reactor was denoted as “small scale Sabatier reactor”. The model assumptions and 
equations were slightly different in each Sabatier reactor, and the differences will be discussed in the “Model 
Assumptions and Equations” section. 


A. Small Scale Sabatier Reactor Description 

The small scale Sabatier reactor consisted of a 304L stainless steel tube that was 1.25 cm OD and 1.1 cm ID. The 
reactor was packed with two different bed lengths of the same Ru/A1203 catalyst pellets as used in the APM Sabatier 
reactor. The bed lengths tested were approximately 4.5 cm and 8.89 cm long, yielding catalyst volumes of 2.3 ml (2.18 
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g) and 6.3 ml (6.67 g) respectively. The schematic of the reactor and different bed lengths is shown in Figure 5, with 
the catalyst bed located in the shaded/crosshatch area of the drawing. The 4.5 cm bed matched closely with the 
diameter of that found and modeled in literature’, which provided additional literature for data validation. Five TCs, 
labeled A — E, were placed on the outside wall as shown in the figure. TC F was placed in the center of the catalyst 
bed and TC G was at the end of the catalyst bed. The reactor was sealed with Swagelok fittings and wrapped with heat 
tape (for pre-heat) and a layer of insulation. 

The reactor was operated by pre-heating the catalyst bed until it reached ~250 °C under a Hz flow. The gas flow is 
from left to right in the schematic. The CO2 flow was started and once the reactor temperature increased, the heating 
tape was turned off. The H2 to CO: ratio was 5:1 in all of the tests, and multiple flow rates were varied at each reactor 
volume, with the reactor pressure at 2.7 atm in all experiments. The reported reactor temperatures and CO2 conversion 
percent were measured after the reactor temperatures had stabilized (approximately 1 to 3 hours of operation). 


22cm 2.2cm 


8.9cm 
45cm 


Figure 5. Left: 2.3 ml of catalyst; Right: 6.3 ml of catalyst. 


IV. Model Assumptions and Equations 


The initial section of the model description describes the assumptions and equations for the APM Sabatier reactor. 
The last section will describe the modifications made for the small scale Sabatier reactor. The Sabatier reactor model 
was first developed using Python software. The model aimed to characterize the Sabatier reactors in terms of tracking 
the temperature and the ratio of the various components in the mixture throughout the reactor (i.e. mole fraction, yj). 
The first iteration of the model was simple, while future versions will add complexity in order to increase accuracy. 
The Python script worked by breaking the reactor domain into incremental pieces and solving equations within each 
segment. Figure 6 shows an example of a single segment along the reactor volume (AV). The current model 
maintained a stoichiometric ratio of 4:1 for the reactant gases H2:COz for validation, with the ability to modify this 
ratio in the future. The initial assumptions made on the system were as follows: 


Plug flow conditions were in the reactor bed (no radial profile). 
Axial diffusion in the reactor bed was neglected. 

The system was at steady state. 

The outer wall of the reactor was adiabatic. 

No side reactions occurred. 

Ideal gas laws applied. 

Uniform porosity existed throughout the catalyst bed. 

The effectiveness factor, n, was 1. 


av 
Figure 6. Reactor domain segment. 


GOS ON: Pe Ga Bo 


All known variables (flow rate, inlet temperature, initial species, etc.) were defined at the beginning of the program. 
The main function of the program is outlined in the flow diagram in Figure 7. This function determined the differential 
equations that define the change in species, reactor temperature, feed temperature, and pressure drop based on the inlet 
conditions for a given domain segment. The integration of these equations across the entire volume determined the 
conversion and temperature profile of the reactor. 

The “Calculate Mass-Weighted Fluid Properties” function used a library of thermophysical properties that 
calculated fluid properties for the mixture of chemical species at the specified reaction temperature based off of the 
mass of each chemical species present. The fluid properties were taken from a NIST database.’ The fluid properties 
included Cp, p, k, and yw as well as an updated average velocity of the fluid based on the density change. Cp was 
calculated using polynomial coefficients.!° The specific heat polynomial coefficients for enthalpies of formation were 
gathered from literature.'! The fluid properties were used to determine the effective overall heat transfer coefficient; 
this determined the amount of heat transferred between the reactor and feed which drove the change in temperature of 
each throughout the volume. A change in pressure was calculated based off of the updated fluid properties as well. 
The species composition changed based on the reaction rate for the given pressure and temperature conditions. 
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Figure 7. Main Python programming function flow diagram. 


A. Kinetics 

Some literature has modeled the rate of the overall Sabatier reaction via the Hougen-Watson rate equation |” as 
well as empirically '°. The empirically derived method was used in this work which modeled the rate of disappearance 
of COz and is shown in Eq. (2). 


r = ko exp(—E,/RT) fc" = Wen] Pasol “Prnel po /(er) (2) 


B. Material Balance 

The steady state mass balance of the reactor “ was calculated using Eq. (3), and conversion of CO2, X, was 
calculated using Eq. (4). The internal effectiveness factor is a measurement of how far the reactants diffuse into the 
pellet before reacting. In essence, it measures how effectively the catalyst is being used. When it is less than 1, it is 
diffusion limited, and when it is near 1, the pellet is reaction limited.!3 In this work, it was assumed that n =1. The 
intrinsic kinetics on the surface of the catalyst pellet have not been studied in great detail for non-isothermal reactions. 
We hope to study the non-isothermal conditions and effect of diffusion limitations by the catalyst in future work. 
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dN; 
aT = yar (3) 
N 
X=1- 602, (4) 
c02,0 


C. Pressure Drop 

For slim packed tubes, wall effects are an important factor when considering the pressure drop across the packed 
bed. An appropriate correlation for a packed bed with cylinders was given by Reichelt.!° The pressure drop equation 
is given in Eq. (5). 


dP MW [190A,,7ug», (1 —€)? Ayu 1 
OF 190A Warglt Gl 8)" Avtowal ai /101325 (5) 
dV Ac D, = ByDy 
where, 
2 D.\2 . 
A, = 1-4. ——— B, = 2.00 (—") + 0.77 
“ 3(D,/D,)(1 — €) us D;, 
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D. Energy Balance 

The amount of heat transferred between the reactor and feed gases determined the temperature profile through the 
reactor domain. The methodology by which it was calculated was the thermal resistance method. An overall heat 
transfer coefficient was calculated based on the average effective radial thermal conductivity of the bed and the 
convective heat transfer along the wall, resulting in Eq. (6)!°. 


Upp ih 6 


1 1 1/Bi+3\ D 

2 ( )— (6) 
The convective heat transfer coefficient is calculated by the Nusselt number correlation as suggested by Martin & 
Nilles'’ shown in Eq (7). 


hD 
Nu= s 


D,\ (k 
= (13 cy 5°) ( “<) + 0.19Re®75pr0-33 (7) 
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The effective radial thermal conductivity is given by Eq. (8) based on correlations by Winterberg.!*!? The radial 
average of Eq. (9) was used in the calculation of Eq. (8). 


1 
Ref = Kpea + are fZ2.-ryk (8) 
2 
Z,-1r' 
eee ,  2£1<Z,—-17' <D,|0.44 + 4exp(—Re 9 
a (0.44 + 4exp(—Re/49)) _ ! tae a 
1, Otherwise 


20,21 


The stagnant bed conductivity based on work by Zehner and Schlunder was calculated from Eq. (10). 


Races: ae ~e[G- 8B | 1 lane (10) 
per ae ale ~EBy! n(s)- @+-25] 
k 1l—e 10/, 
Wh = d = 
ere é Ea an B= 2.5( ; ) 


The constant calculated from Eq. (11) was used to calculate the heat transfer along the interface between the reactor 
and internal feed pipe. It is a ratio of the surface area of the feed pipe to the reactor cross-sectional area. 


Spt 62.19 
~ OD2—ID2 (11) 


The energy balance for the reactor for a single reaction and only one heat transfer boundary, is shown in Eq. (12). 
When integrated across the volume, this equation determined the temperature of the reactor. 


dT _—AHgr + Veep B(Ty — T) 


one (12) 
dV di NiCp; 


Because there was no reaction within the internal feed tube, the energy balance in the feed tube only included heat 
transferred between the reactor and feed gas, as described by Eq. (13). 


aT, U Ty —T 
7 _ Vers B(Tr -T) Where 
av Di Nip Pir 


Trlv=vp = Tro (13) 
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In order to find steady state conditions, a “shooting reactor out reactor in 


strategy” was used. An initial guess for the reactor inlet feed in + feed out 
temperature was provided to the model, and integration —— = 
over the reactor volume yielded an inlet feed temperature V=Vp V=0 

value for the specified reactor inlet temperature. This Trlv=v_ = Tro Tly-0 = Trlv=o + Taips 
value was compared to the known inlet feed temperature 150 Steady State Solution 


and iterated until there was a match. A sample output from ~~ Tp 
this iteration of reactor conditions is shown in Figure 8. 2007) = Tilv=vn | 
The intersection points indicated a steady-state condition. 
There is a caveat to the boundary condition at the reactor 
inlet/feed outlet. Based on experimental results, there is of 
some heat leak through the preheat loop as a result of 


imperfect insulation. While the amount of this heat leak is an | 
difficult to quantify, the model attempts to incorporate this 100 | 
based on a specified temperature difference (Tair) between 

the feed outlet and the reactor inlet. -150 | 


E. Small Scale Sabatier Reactor Equations 


The same kinetic and material balance equations (Eq. _250 | | | | | 
Se : 0 100 200 300 400 500 600 
(2) through (4)) were utilized for the small scale Sabatier Thv-0 FO 
reactor along with the same pressure drop equation (Eq. Figure 8. Top: Reactor domain boundary conditions; Bottom: 
(5)). While a major difference between the two reactors Steady-state determination for inlet reactor temperature. 


was the lack of external cooling and internal feed tube, the 

same energy balance equation (Eq. (11)) was also used. The results section will discuss why this term was included 
in the energy balance. The model can predict the results of adiabatic operation by setting the effective heat transfer 
coefficient to zero. 


V. Results and Discussion 


A. APM Sabatier Reactor Experimental and Modeling Results 

Long duration (7 hour) runs were performed on the reactor to determine steady state reaction temperatures during 
nominal operation of CH, production. The results of the thermocouple temperatures at steady state for nominal 
conditions (with Hz recycle loop on) and for single pass conversion (no H2 recycle loop) are displayed in 

Table 1 and Figure 9. At nominal conditions and steady state, >99% conversion of CO2 to CH4 was performed and 
validated via GC, along with 32 grams per hour of H2O production. TC14-C was the highest temperature of the preheat 
tube, at 327 °C for nominal conditions, and TC 8, near the reactor inlet, was the hottest spot of the entire system. Due 
to thermocouple spacing, it was possible that there were regions at higher or lower temperature than recorded by TC 
locations. The catalyst bed reactor inlet, and catalyst bed reactor outlet were 260 °C different in temperature during 
for nominal conditions. By the time the gases left the reactor outlet, they were cooled down to the point where the 
Sabatier reaction was not expected to have much activity and were also interfacing with the coolest part of the preheat 
tube (87 °C at TC 14-A), which also decreased the temperature. The reaction should have the highest conversion near 
the catalyst bed inlet and then subside as it travels out of the reactor bed and leaves the system. Because of the lower 
temperature profile (<300 °C), the reactor may be oversized for nominal flow rate conditions, but there is a threshold 
of higher throughput flow rates, as observed with catalyst damage occurring at 1250 sccm CO2 and 5000 sccm H2. 


Table 1. Temperature averages at steady state for a nominal run, displayed in °C 
HL/CO; Recycle | TC14-A | TC14-B) TC14-C TCI5 TC7 
TC8 TC6 TC9 TC10 
9 
flow rate, loopon? (Preheat (Preheat | (Preheat (loop (loop (inlet) | (mid.) | (mid. (out) 


Single Pass 
Co, 


sccm Y/N 42 cm) 27 cm) 20 cm) inlet) exit) conversion, % 

750/3000 Y 87 192 327 142 132 425 264 254 165 - 

250/1000 N - - 431 59 54 448 183 174 129 95 
350/1400 N - - 469 71 57 502 205 186 117 92 
450/1800 N - - 505 93 70 541 252 227 121 92 
550/2200 N - - 527 113 82 564 297 270 129 91 
650/2600 N - - 539 132 95 575 335 307 135 88 
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For single pass Joutuet 
conversion, the flow rate 
was gradually increased 
from a H2/COz flow rate of Inlet 
750/3000 sccm to 650/2600 ~~ 
sccm. As the flow rate was 
increased, the conversion of 
CO> in the system declined. 

The addition of the recycle — 
stream allows more Figure 9. Temperature profile through APM Sabatier reactor at steady state conditions. 
complete conversion of the 

CO>. This could be due to the fact that a small additional stoichiometric amount of H2 is available, and also additional 
gas from the cool recycle stream helps control thermal stability in the catalyst bed. 

The results of the single pass Python model for temperature profiles across the normalized reactor volume are 
shown in the left of Figure 10. The input CO/H2 flow rates of 650/2600 sccm were used for these images. The feed 
tube temperature inside the reactor was manually input from experimental data. The experimental results of the catalyst 
bed thermal profile (blue dashed line), labeled T, and were much lower at the inlet than the model predictions (blue 
solid line). At >600 °C, the catalyst would become damaged and the reaction would shift its hot spot location towards 
the middle and outlet regions of the reactor catalyst bed to achieve conversion. The model predicted a decrease in 
temperature across the volume of the reactor, which was observed experimentally as well. Ts; was the reactor feed 
temperature, which at this case experimental values were inserted into the model. Some adjustment is needed in the 
thermal prediction of the model for the preheat loop and inlet conditions, as the area of heat loss in this region is not 
entirely accounted for in the code. It is expected that as pellet considerations are also modeled, the predictions will 
become more 
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displayed as mole fractions during nominal conditions of the APM 
Sabatier reactor across the normalized volume is show in Figure 
11. The trends also match empirical data for a 4:1 H2:CO2 fed 
reaction and follow what would be expected of the single pass 
conversion of 88%. 
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catalyst surface. Current equipment is limited in how many TCs can be added, and thermal hot spots may be undetected 
at this time. The model will help determine these expected hot spot regions, and TCs can be placed appropriately. 


B. Small Scale Sabatier Reactor Experimental and Modeling Results 

The experimental results of the small scale single pass Sabatier reactor temperature profile is shown in Table 2. 
98% conversion was achieved at the lowest flow rate of the 2.3 ml catalyst bed volume, which had the highest 
temperature of 376 °C at the middle of the catalyst bed, while maintaining the lowest maximum temperature of all 2.3 
ml runs. Total flow rates less than 240 sccm were not evaluated, because at this low flow, the GC did not sample 
properly. As flow rate increased, the temperature in the catalyst bed increased, but conversion decreased. For the 
6.3 ml catalyst bed, the same trend of decreased conversion at increased flow rate was observed. Since increasing the 
volume of the catalyst bed in this reactor resulted in a longer bed of the same diameter, the lowest flow rate evaluated 
was selected such that it was slightly lower than the maximum flow rate in the smaller volume reactor. Since 376 °C 
is roughly the high temperature needed for >98% conversion to CH4, we can discuss possible scenarios for Mars 
Sabatier reactor development. If several small reactors can be placed in series and run at lower flows, the thermal 
integrity of the catalyst is protected, and >98% conversion is ideal. If one or two reactors are needed at high throughput 
of gas, the system may perform at lower conversion of COz, and also run the risk of thermal damage to the catalyst 
pellets (>600 °C). 


Table 2. Experimental data from small scale Sabatier reactor. 


Renee tee - oe Thermocouple, °C Conversion 
A Cc E F G 
ml sccm surfs B surfa D ; : atalys x 
Hey. mide) (surface outlet) | og mide) |_bed outlet 
2:3 40/200 225 295 314 295 258 376 293 98 
2.3 80/400 106 236 289 259 175 417 259 96 
2.3 120/600 93 282 366 333 219 523 337 97 
2.3 160/800 85 298 416 395 267 570 407 93 
6.3 150/750 347 391 362 295 249 399 252 98 
6.3 165/825 350 408 383 314 264 424 271 99 
6.3 225/1125 364 455 454 398 345 501 355 97 
6.3 275/1375 361 472 484 440 390 531 405 91 
6.3 300/1500 355 478 494 456 410 542 427 85 
6.3 350/1750 336 486 511 481 437 558 460 91 
6.3 375/1825 329 485 516 489 449 564 471 82 
6.3 425/2125 253 467 530 512 475 579 500 80 


The Python model results of the small scale single pass Sabatier reactor is shown in Figure 12 and Figure 13.The 
images represent the larger bed volume (6.3 ml) and a CO2/H;2 flow rate of 425/2125 sccm. If the operation is assumed 
adiabatic, the results indicate a very poor performance of the reactor as shown by the temperature and conversion 
profiles in the left side of Figure 12. By implementing the wall temperature profile as obtained from the experiment, 
a much more accurate set of results are predicted by the model as can be seen in the plots to the right in Figure 10 and 
Figure 13. With the adiabatic case, the entire reaction stops when it reaches the maximum temperature. It is apparent 
that some cooling is necessary to maximize the conversion within the reactor and that the experimental reactor is not 
as well insulated as assumed. The model predicts the hot spot is located approximately 25% into the reactor. In 
addition, the model for the non-adiabatic case shows that highest temperature in the bed exceeds 650 °C, which is 
thought to be high enough to damage the catalyst. This shows how the model can be used to determine that a particular 
reaction condition is not appropriate for the reactor. 

The catalyst is above 600 °C for only about 1 cm of the reactor length. There can be only a finite number of 
thermocouples in the catalyst bed, and it is possible that the hot spot would be missed as occurred in this case. 


C. APM Sabatier Reactor Catalyst Characterization Results 


I. XRD Results 

The XRD scan of the entire Ru/Al,O3 pellet from an unused pellet, a damaged pellet from 2015, and the current 
2017 Sabatier reactor with a pellet from the reactor inlet, middle and outlet pellets is shown in Figure 14. The low 
intensity and broad peaks reveal that the material is amorphous. Since the diffraction peaks remain unchanged after 
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Figure 12. Python model results for the small scale Sabatier reactor. Left: Adiabatic operating condition; Right: Imposed wall 
temperature profile based on experiment. 
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Figure 13. Species mole fraction profiles for the small scale Sabatier reactor. Left: Adiabatic operating condition; Right: Imposed 
wall temperature profile based on experiment. 


and during use at elevated temperatures, no significant change in the crystallinity occurred pre- or post-reaction. The 
pellets show reflections corresponding to AlzO3 (A) and RuO2(e). a-Al2O3 and y-Al203 are the most common phases 
of Al2,O3. a-Al203 is polycrystalline and typically used in ceramic applications, is the most stable, and calcined at high 
temperatures (>1000 °C), whereas y-Al2O3 has high surface area, processed at lower calcination temperatures (~500 
°C), and is typically used in catalyst applications. It is often observed in literature that that the y-Al2O3 transforms into 
the a-Al,O3 phase at elevated temperatures.”” It was confirmed by the catalyst vendor that the Sabatier pellets 
purchased for this were supported on y-Al203. 
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Figure 14. XRD Spectrums of Ru/Al,0; pellets from APM project, unused, used in 2015 (damaged), and used 2017 (not damaged). 


The AloO3 phases were observed at 37.8°, 39.5°, 45.5°, and 67° 20. The 37.8° and 67° 20 were stronger in o-Al203, 
corresponding to the 110 and 214 reflections, respectively. Although faint, the value at 28.4° was most prevalent in 
the unused, catalyst bed middle, and catalyst bed outlet pellets, which corresponded to the RuOz 110 reflections. The 
unused catalyst has not been heat treated yet, and the catalyst bed outlet sustained the lowest temperatures during 
operation, so it was not participating in the high exothermic reaction conditions and likely not undergoing thermal 
shock like the inlet and middle catalyst bed were. At 39.5° and 67°, the 20 values were also characteristic of the RuO2 
200 and 301 reflections, respectively, but also overlap with the Al2O3 reflections at these locations. RuO2 species was 
likely dispersed well on the Al2O3 support surface prior to high temperature (>200 °C) exposure. Confirmation of 
metal sintering or particle aggregation was challenging to determine via XRD since the pellets were ground in the 
aqueous phase during preparation for analysis. Even if sintering did occur at high temperature reaction conditions, the 
size of crystallite may not have grown large enough for XRD detection.” 

The XRD scan of the Ru pellet powder and Al203 pellet powders from the current 2017 in-use Sabatier reactor 
system (used but not damaged) is shown in Figure 15, along with reference spectrums of the y and a Al2O3 phases, 
Ru, and RuO». The main difference observed from these scans was the fact that the cooler portions of the catalyst bed 
(the outlet and middle) had reflections corresponding to RuOz, while the hotter reactor inlet did not. 

It was challenging to completely separate the powdered Ru coating from the Al203, and therefore some overlap 
was present when attempting to isolate the two materials. The reference reflections also overlap with the experimental 
reflections for the Ru, RuO2, and Al2O3 patterns, especially for broad peak results. The reference o-Al2O3 matched 
closest with the support structure, but some of the reference y-Al2O3 peaks were also present. 


2. XPS Results 

The XPS Ru3D scan for the Ru/Al,O3 pellet of an unused pellet, damaged pellet from 2015, and inlet, middle and 
outlet pellet of the current 2017 APM reactor (used but not damaged) are shown overlayed in the top of Figure 16. 
The left image (top of Figure 16) was at the surface (etch level = 0) and the right was at approximately 15.75 nm 
below the first oxide layer (etch level = 9). The chemical states of the Ru3ds5,2 spin orbit is 280.2 eV for Ru Metal and 
280.7 for RuO>.”4 On the surface (etch level = 0), the inlet, middle, and outlet of the 2017 APM reactor pellets were 
similar with binding energies of the right at ~280.2 and the left at ~284.6 eV. The unused pellet had the farthest shift 
to the left, and the used/damaged pellet had the farthest shift to the right. The same behavior was observed after etching 
away some of the surface oxide (etch level = 9), but the used pellet had similar binding energy to the 2017 pellets. 
Since y-Al203 has higher surface area, the pores in the catalyst support may have collapsed during the phase change 
to a, encapsulating the Ru, which could be why it was not present in the hotter regions of the reactor, and the overall 
surface area of RuOz would decrease. This damage may have contributed to the catalyst inability to become 
successfully reduced for CH selectivity.”* 
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The XPS Al2p scan for the Ru/AlO3 
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VI. Future Work —_ 


Computational Fluid Dynamics (CFD) Figure 15. XRD Spectrums of Ru and AlI,O; isolated powders from the pellets 
simulations are currently being conducted from current APM Sabatier reactor. 
using STAR-CCM+ as a separate way to 
verify modeling results. The CFD approach 
allows the ability to incorporate 3D flow 
effects into the model. Because the ratio 
between the tube diameter to pellet size is low, 
it is important to include the non-uniformity 
of the flow and the wall bypassing effects that 
it may cause. The validity of certain 
simplifying assumptions in the Python model 
can be determined by comparing the results of 
the CFD to the Python script. Once the Python 
script is at an acceptable accuracy, it can be 
used to determine other suitable reactor 
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approach. Figure 7 shows a close-up of the Figure 16. Top: Ru3d Overlay. Left: Etch Level 0, Right: Etch Level 9; 

reactor domain mesh with the fluid in blue and Bottom: Al2p overlay. Left: Etch level 0, Right: Etch Level 9. 


the solid catalyst pellets in black. 

There are also plans to expand the capability of the Python model to 
take into consideration the energy and mass balance across the catalyst 
pellets and effects of non-isothermal particle considerations. Since the 
reaction is exothermic, the catalyst pellet internal temperature may rise. We 
must keep track of two phases being affected by this reaction: the fluid- 
phase and the diffusion process on the solid-phase catalyst pellets. These 
phases are simultaneously exchanging mass and energy and will be 
modeled for their relationship to one another. The energy balance equations 
will be coupled into the enthalpy change on the fluid, and the enthalpy 
change on the particle. The Weisz-Prater criterion can be used to measure 
the rate of reaction to determine if there are any internal diffusion 


pane : pest : Figure 17. CFD reactor domain mesh 
limitations of the reaction. The Mears’ Criterion can be used to observe if with fluid (blue) and solid catalyst 


external diffusion is limiting.!* (black). 
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VII. Conclusion 


A thermal and kinetic model of the APM Sabatier reactor and small scale Sabatier reactor was created using 
Python. The model was developed after a catalyst performance failure occurred in the APM Sabatier reactor in 2015 
during higher than nominal flow rate testing. Since high flow rates will damage the catalyst, the model will help 
optimize the throughput of the catalyst bed, and be validated experimentally to build trust in the model. The model 
will continue to evolve including considering transport by convection in the fluid, diffusion inside the catalyst pellets, 
reaction kinetics, and thermal effects of the gas-phase and solid-phase aspects of the reactor system. The reaction 
diffusion investigation of the catalyst will be coupled to the mass and energy balances of the fluid, as well as the 
pressure drop considerations across the reactor. The continued development of this modeling will allow us to optimize 
systems that can predict and determine acceptable flow rates and temperature profiles in order to avoiding damaging 
the catalyst and troubleshoot system tolerances. It will also help with scaling larger systems or systems in series for 
future design of a Mars ISRU propellant production system. 

Ru/AlO3 catalyst characterization analysis resulted in subtle changes that were observed in both XRD and XPS. 
The XRD crystal structure on the pellets that encountered higher reaction temperature (>500 °C) were missing RuO2z 
characteristic reflections, and the cooler portions of the catalyst bed (reactor middle and outlet) had reflections 
corresponding to RuO2, while the reactor inlet did not. The XPS observations revealed that a stronger oxide layer was 
present on ‘healthy’ or unused catalyst as well, but not present on the damaged catalyst. If the phase change from y- 
Al203 to o caused pores to collapse and obstruct the surface layer of Ru, this may have been why performance for 
Sabatier selectivity and reaction temperatures declined. Overall, catalyst exposure to reaction temperatures greater 
than 450 °C caused the thermal shock on the catalyst surface, losing the ability to selectively form CH,. 
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